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ABSTRACT 



We point out that in compact U(l) lattice gauge theory on a 
finite volume, some topological excitations (monopole loops and 
Dirac sheets winding through the lattice) effectively divide 
the configuration space in separate sectors. Local Monte Carlo 
algorithms have problems in moving from one such sector to 
another and may thus produce misleading results. We have ob- 
served this phenomenon in two independent programs, one for the 
Z(128) subgroup and Villain's action near the phase transition, 
the other for Z(32) and Wilson's action in the cold phase. 
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We observed very long- living metastable configurations in MC-simulations 
of d=4 U(l) gauge theory. As reason for this behaviour we identified objects 
specific for lattices with finite size topology. These are on one hand mono- 
pole current loops (world lines) winding through the lattice and closed 
only due to the periodic boundary conditions. On the other hand we found 
Dirac sheets, the world sheets of Dirac strings in three dimensions, without 
monopole loops as boundary and closed due to periodic b.c. 

Both objects are responsible for long living action gaps: the monopole 
loops at the PT, the Dirac sheets throughout the cold phase. They may 
obscure the nature of the PT and other properties of the system. This 
motivates the presentation of our results. 

The importance of monopoles for the dynamics of the PT has been discussed 
in the literature [1,2,3,4,5]. They have been observed in MC simulations of 
U(l) lattice gauge theory and found to correlate strongly with the free 
energy at the PT[4]. In d=3 monopoles may be regarded as endpoints of 
Dirac strings of magnetic flux. They are identified, using Gauss' law, by 
measuring the total flux leaving a closed surface like e.g. the boundary 
of a cube. In d=4 one proceeds in an analogous way and we follow Ref.2 
for the definition. We introduce plaquette flux variables 9 pa ( x ) as the 
sum over the link angles 0 € [ -HIT ] 
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and the physical flux 
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Fig. 1 shows the accumulated total length L of monopole loops 

L(x) = I nil, < 8 
£> x * 

where n is the number of occurences of loops with length %. We also 
indicated the steps due to the contributions from the loops, that close 
because of periodicity only. Subtracting such contributions from the curve 
for configuration (a) one obtains a curve very similar to that of configura- 
tion (b) (cf. Fig. lb). 

Due to our periodic and untwisted b.c. there can be no single periodically 
closed monopole loop, it has to have at least one partner in the opposite 
direction. The total number of boundary crossings, taking all periodically 
closed monopole loops together, is therefore zero. 

Although the fluctuations at the PT are of long range the loops winding 
around the lattice have difficulties to decay. Their only chance is to 
cooperate with another such loop in order to undo the surplus boundary 
crossings. For local MC-algorithms this is very unprobable and thus it 
may take a very long time - hence the apparent stability. 

In order to investigate the decay properties the last configuration (a) was 
subjected to further MC-iterations at 3 = 1. Within the first 50 iterations 
the total length L(0) decreases from 14880 down to 2850. This is mainly 
due to the vanishing of ordinary closed monopole loops. For 900 subsequent 
iterations one observes essentially just one pair of (periodically closed) 
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monopole current loops, which do not succeed to annihilate. The last 
configuration considered had still total length L(0)=456 and only this 
monopole pair left. Had we cooled the system "adiabatically" the monopoles 
might have had a better chance to drop out. 



What happens to the Dirac sheet when a monopole loop closes and dis- 
appears? On an infinite lattice without the periodically closed loops the 
sheet must disappear as well. On a finite lattice, however, there are 
two possibilities. In Fig.2 we illustrate the situation on a 2-dimensional 
plane of the dual lattice A (for linear size 8). Situation (b) is what 
happens also on an infinite lattice. In case (a) one is left over with a 
trapped magnetic flux, a Dirac sheet covering a whole surface closed 
due to the periodic b.c. Details of the shape of the sheet are not gauge 
invariant but the fact of its existence is. 



In independent MC-simulations for Z(32) on lattices of various size and 
for Wilson's action, 



^ W yv _ exp(-0(l-CGse uv )). (9) 



we observed metastable states of enormous lifetime. They occurred 
frequently at relatively large 3 -values (in the cold phase), when the start 

configuration was disordered (hot). The smallest observed action gap varied 

4 4 
with the lattice size and was e.g. 0.0127(4) for 4 , 0.0026(2) for 6 , and 

0.0007(1) for 8^ lattices for 3 =3.5. The 3 -dependence is weak as may be 

seen in Tab.l. For that gap no decay from the excited state to the ground 
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state was observed in over 12000 iterations for $>1.1. 

The gap was noticed originally in thermal cycles of a gauge-Higgs system 
[8 ] for fixed 3 and varied gauge-Higgs coupling k . The plaquette remained 
excited in the upward branch (increasing k ) of the hysteresis until far 
above the Higgs transition. On the downward branch the plaquette 
variable usually selected the ground state and remained there in further 
cycles. As this phenomenon occurred in the gauge variable it was not 
related to a similar effect observed in Ref. [ 9 ]. 

In order to identify the reason for the metastability we performed an 
iterative gauge transformation procedure for an excited configuration as 
well as for a ground state configuration. At each site a gauge transforma- 
tion 



6 T ^fl' = 8 ¥4 ± xJmod 2n) * [-11,11] 

X,±y X>±|j x i±\X x 



(10) 



is chosen such as to minimize 



I I9' x>u l <»> 



In ambiguous situations the value of is chosen randomly. Notice that 
a gauge transformation even with small X v can change the value of 
n p of a plaquette by flipping some link angle e.g. from a value close 
to - II to a value close to TL Such a gauge transformation is applied to 
all sites of the lattice which defines one "gauge iteration". The procedure 
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tends to decrease the number of plaquettes with n^O (plaquettes with 
negative n p may cancel others with positive n p ). After sufficiently many 
gauge iterations a minimal number of such plaquettes is obtained. Fig.3 
shows the yz-plaquettes with n p £0 (i.e. possibly contributing to a Dirac 
sheet) as found for an excited configuration (a) and a ground state con- 
figuration (b) after zero, one and six gauge iterations. Whereas the 
distribution of plaquettes with n p ^0 before the gauge iterations is quite 
similar for the excited and the ground state configuration, one finds a 
remarkable difference after several gauge iterations. 

Note that the total number of plaquettes with n^AO for each xt-plane is 
always odd for the excited state but even for the ground state. The 
excitation is identified with a Dirac sheet covering a whole torus on the 
dual lattice. 



Let us estimate the resulting action gap for the case of Wilson's action 

in the weak coupling situation ( 3^°°). For a Dirac sheet in the xt-plane 

of the dual lattice the minimal surplus flux through every zy-plane of 

the original lattice is 2 II due to one plaquette with n p A0. Eventually it 

2 

will be distributed over all N plaquettes of each zy-plane (for lattice 

4 2 
size N ). Each of these N plaquettes contributes in average 

1 - cos(2lI/N 2 ) ffc 2n 2 /N 4 (12) 
leading to a minimal gap of in the plaquette energy averaged 
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FIGURE CAPTIONS 



Fig.l (a) The integrated total length of monopole loops L(x) for a 
configuration (denoted by a) with higher and a configuration (denoted 
by b) with lower average plaquette energy E at 3=0.643; the arrows 
indicate the steps due to contributions from periodically closed 
monopole current loops. 

(b) Here we plot again the same curves as in (a) on a different 
scale after subtracting off the contributions from periodically closed 
loops from curve a. 

Fig.2 We illustrate the possible history of a pair of periodically closed 
monopole loops depicted on the dual lattice. 

(a) Here a Dirac sheet (shaded area) is left over after the disappear- 
ance of the loops, and 

(b) no Dirac sheet is left. 

Fig.3 We compare an excited (a) and a ground state (b) configuration 
(obtained at 3 =3.5 from hot and cold starts after several hundred 
MC-iterations); the x, y, z, and t coordinates may be identified 
as indicated, the squares denote zy-plaquettes without (open square) 
and with (full square) 2 II -flux running through. The first column 
gives the configurations before gauge iterations, the second and 
third column after one and six gauge iterations, as discussed in 
the text. 
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